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Following  the  theory  of  error  propagation,  we  developed  analytical  functions  to  illustrate  and  evaluate 
the  uncertainties  of  inherent  optical  properties  (lOPs)  derived  by  the  quasi-analytical  algorithm  (QAA). 
In  particular,  we  evaluated  the  effects  of  uncertainties  of  these  optical  parameters  on  the  inverted  lOPs: 
the  absorption  coefficient  at  the  reference  wavelength,  the  extrapolation  of  particle  backscattering  coef¬ 
ficient,  and  the  spectral  ratios  of  absorption  coefficients  of  phytoplankton  and  detritus/ge  lbs  toff,  respec¬ 
tively.  With  a  systematically  simulated  data  set  (46,200  points),  we  found  that  the  relative  uncertainty  of 
QAA -derived  total  absorption  coefTicients  in  the  blue-green  wavelengths  is  generally  within  ±10%  for 
oceanic  waters.  The  results  of  this  study  not  only  establish  theoretical  bases  to  evaluate  and  understand 
the  effects  of  the  various  variables  on  lOPs  derived  from  remote-sensing  reflectance,  but  also  lay  the 
groundwork  to  analytically  estimate  uncertainties  of  these  lOPs  for  each  pixel.  These  are  required 
and  important  steps  for  the  generation  of  quality  maps  of  lOP  products  derived  from  satellite  ocean  color 
remote  sensing.  ©  20 10  Optical  Society  of  America 
OCIS  codes:  010.4450,  280.4991. 


1.  Introduction 

To  describe,  understand,  and  forecast  the  interac¬ 
tions  and  changes  associated  with  environments  and 
climate,  it  is  critical  to  have  reliable  long-term  obser¬ 
vations  of  geophysical  properties  for  the  atmosphere 
as  well  as  for  the  land  and  oceans.  For  such  a  daunt¬ 
ing  requirement,  measurements  by  satellite  remote 
sensors  are  indispensable.  To  meet  this  requirement, 
systematic  observations  of  the  biogeochemical  prop¬ 
erties  of  global  oceans  through  the  measurement  of 
ocean  color  (ocean  color  radiometry — OCR)  from 
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space  have  become  an  important  component  of  the 
Earth  Observing  System  [1],  and  various  ocean  color 
missions,  such  as  the  CZCS  (Coastal  Zone  Color 
Scanner,  1978-1986),  SeaWiFS  (Sea-viewing  Wide 
Field-of-view  Sensor,  1997-present),  MODIS  (Moder¬ 
ate  Resolution  Imaging  Spectroradiometer,  1999- 
present),  and  MERIS  (Medium  Resolution  Imaging 
Spectrometer,  2002-present)  instruments  that  have 
been  launched  since  the  late  1970s  [2]. 

The  derived  products  from  OCR  generally  include 
inherent  optical  properties  (lOPs;  e.g.,  marine  ab¬ 
sorption  and  scattering  coefficients),  concentrations 
of  chlorophyll-a  and  suspended  matters,  and  water 
clarity  [8].  These  products,  however,  inherently  con¬ 
tain  some  degree  of  uncertainty  because  of  imperfect 
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sensor  engineering  and  retrieval  algorithms  [1-7]. 
Conventionally,  such  uncertainties,  sometimes  called 
“errors,”  are  evaluated  statistically  by  comparing  re¬ 
trieved  values  from  remote  sensing  with  those  from 
water  samples  of  concurrent  measurements,  and  an 
averaged  uncertainty  for  the  entire  data  set  [8],  or 
subgroups  [0],  is  usually  obtained.  This  quantity  pro¬ 
vides  a  general  picture  about  the  consistency  be¬ 
tween  retrieved  and  measured  properties.  By  the 
nature  of  the  evaluation  method,  however,  this  quan¬ 
tity,  especially  the  averaged  “error”  derived  from  a 
data  set  covering  a  wide  dynamic  range,  does  not  re¬ 
present  the  uncertainty  for  a  particular  pixel  in  sa¬ 
tellite  imagery.  Also  note  that  this  quantity  does  not 
necessarily  represent  algorithm  or  satellite  retrieval 
error,  as  errors  associated  with  various  in  situ  or  lab 
measurements  are  not  removed  in  the  comparison. 

At  present,  there  is  no  quantitative  measure  able 
to  describe  the  quality  (or  uncertainty)  of  satellite 
products  pixelwise,  although  eight  quality  levels  of 
chlorophyll-a  concentration  could  be  obtained  based 
on  the  class-based  approach  recently  developed  by 
Moore  et  al.  [9].  Consequently,  this  limits  the  uncer¬ 
tainty  evaluation  of  next  level  products,  such  as  pri¬ 
mary  production  ( 1 0, 1  1  ]  or  water  quality  [12-14]  and 
prevents  a  full  account  of  errors  or  uncertainties  in 
closure  analyses  [15]  or  trend  evaluations  [lb]. 

Note  that  satellite  ocean  color  remote  sensing  is  a 
complex  system  with  many  variables  and  steps  in¬ 
volved  before  bio-optical  products  are  finally  gener¬ 
ated.  Major  variables/steps  include,  in  a  sequential 
order,  sensor  calibration  [17],  atmospheric  correction 
[IS],  and  inversion  of  optical  and  biogeochemical 
properties  [s,U)]  from  the  spectrum  of  water-leaving 
radiance  (or  remote  sensing  reflectance).  Each  vari¬ 
able/step  has  its  own  sources  and  levels  of  uncertain¬ 
ties;  thus  it  is  required  to  understand  the  sources  of 
uncertainties  and  quantify  the  quality  of  the  pro¬ 
ducts  of  each  step  [6]. 

For  the  bio-optical  inversion,  Wang  et  al.  [20]  pro¬ 
posed  a  scheme  to  obtain  ensembles  of  solutions  for 
each  spectral  remote  sensing  reflectance  spectrum 
(r^s,  in  units  Sr“^,  representing  subsurface  re- 
mote-sensing  reflectance,  which  can  be  derived  from 
above-surface  remote-sensing  reflectance,  after 
correcting  the  air-sea  surface  effects  [21])  via  a  lin¬ 
ear  matrix  inversion  method  [22].  Depending  on  the 
setup  of  the  ranges  and  steps  of  the  spectral  shapes  of 
the  major  components,  such  an  approach  could  result 
in  more  than  1000  combinations  of  inversions  to  be 
carried  out  for  each  spectrum  [20].  This  computa¬ 
tional  requirement  makes  it  rather  difficult  to  be  in¬ 
corporated  into  operational  processing  where  large 
quantities  of  satellite  images  have  to  be  processed. 

The  quasi-analytical  algorithm  (QAA)  [12,25],  on 
the  other  hand,  derives  lOPs  from  r^  in  a  stepwise 
fashion  with  every  variable  also  explicitly  described. 
This  characteristic  makes  it  easy  to  propagate  and 
evaluate  the  effects  of  uncertainties  of  each  variable 
based  on  the  theory  of  error  propagation  [24].  To  high¬ 
light  this  feature  and  to  evaluate  the  uncertainties 


introduced  by  QAA,  we  developed  analytical  expres¬ 
sions  to  describe  the  uncertainties  of  optical  variables 
and  their  propagations  in  the  QAA  process  (see  Fig.  1 
for  a  schematic  flowchart  of  the  steps  and  the  vari¬ 
ables  whose  uncertainties  are  evaluated  and  propa¬ 
gated).  In  particular  we  evaluated  the  uncertainty 
and  propagation  of  the  following:  the  absorption  coef¬ 
ficient  at  the  reference  wavelength,  the  extrapolation 
of  particle  backscattering  coefficient  from  one  wave¬ 
length  to  another  (represented  by  a  power  coefficient), 
and  the  spectral  ratios  of  absorption  coefficients  of 
phytoplankton  and  detritus/gelbstoff,  respectively. 
The  results  of  this  study  provide  important  compo¬ 
nents  that  are  required  for  the  generation  of  quality 
maps  of  lOPs  products. 

2.  Overview  of  Quasi-Analytical  Algorithm  and 
Uncertainty  Propagation 

A.  Quasi-Analytical  Algorithm  Overview 

QAA  is  an  inversion  algorithm  [22]  that  derives 
wateFs  lOPs  from  a  spectral  r^,  as  lOPs  provide  the 
key  connection  between  in-water  constituents  and 
spectral  radiance  from  a  water  body  [25,2()].  The  fol¬ 
lowing  equations  provide  a  concise  summary  of  QAA 
steps  and  data  flow. 

The  quantity  u  is  defined  as 


u(/l) 


bbi^) 

a{X)  + 


(1) 


(±Ac),^  (±Ac)) 

Fig.  1.  Schematic  chart  to  show  variables  and  steps  (S1-S7)  in- 
volved  in  the  QAA  procedure,  redrawn  from  Lee  et  al.  [Jil).  Vari¬ 
ables  with  uncertainties  (Ux-U^)  discussed  in  this  study  are 
highlighted  with  gray,  while  all  others  assumed  error  free. 
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Here  65  and  a  represent  the  total  backscattering 
and  absorption  coefficients  (units  of  m“^),  respec¬ 
tively.  The  symbol  ^  represents  wavelength  (units  of 
nm).  In  the  current  QAA  procedure  and  other  semi- 
analytical  algorithms  [22,‘J7,28],  u  can  be  derived 
from  measured  through  a  model  [29,20]: 


The  absorption  coefficient  at  2  is  then  calculated: 

a(A)=A(A)(66^(A)+6,^,(A)).  (9) 

When  a(2i)  and  0(^2)  are  known,  a<^(22)  and  aph{^2) 
are  further  calculated: 


=  {^0 +^1  (2) 

Here  go  model  constants,  and  their  va¬ 

lues  vary  with  Sun-sensor  angular  geometry  [21], 
although  they  are  considered  wavelength  and  water- 
property  independent.  Equation  (2)  is  an  approxima¬ 
tion  [22]  from  the  radiative  transfer  equation,  so  the 
expression  is  not  error  free.  However,  because  the  er¬ 
ror  in  Eq.  (2)  is  significantly  less  than  10%  [29]  and 
this  study  focuses  on  the  uncertainties  specifically 
introduced  by  the  QAA  procedure,  the  minor  uncer¬ 
tainty  introduced  by  Eq.  (2)  is  omitted  here. 

For  convenience  in  algebraic  derivations,  quanti¬ 
ties  A  and  B  are  defined,  respectively,  as 


A{A)  = 


1  -  u{A) 
m(/1) 


(3) 


(4) 


When  the  absorption  coefficient  at  a  reference 
wavelength  (Aq  —  555,  551,  or  560  nm  for  SeaWiFS, 
MODIS,  and  MERIS,  respectively),  a(>io)>  is  estimat¬ 
ed  through  [1‘)]  (also  see  httj)://\\  \v\v.ioccg.()rg/g roups/ 
Softwaix'  ()(^A/(^AA  vr).[)df  for  recent  updates) 


/,  ^  \a{h)-^a(h)]-[aw{h)-^aw{h)\ 

adg{h)  = - - - - ’ 

g  ^  qUi)]  -  -  Q«,(Ai)| 

(10) 

Here  the  parameters  f  and  ^  represent  the  ratio  of 
aph{h)/aph{h)  and  adgiAi)/adg{^2)y  respectively.  In 
QAA,  the  value  of  f  is  estimated  empirically  from 
rrs(2),  and  f  is  estimated  with  exp [5(^2“'^! )]•  Here 
S  (units  of  nm’'^)  is  the  spectral  slope  of  the  combined 
absorption  coefficient  of  gelbstoff  and  detritus  [28]. 
The  wavelengths  Ai  and  A2  are  normally  selected  as 
412  nm  and  443  nm,  respectively,  to  reflect  spectral 
configurations  of  satellite  sensors  (e.g.,  SeaWiFS/ 
MODIS).  The  subscripts  u;,  p/i,  rfg,  and  p  stand  for 
water  molecules  (including  salt  effects),  phytoplank¬ 
ton,  detritus/gelbstoff,  and  particles,  respectively. 

B.  Uncertainties  and  Their  Propagations 

In  general,  for  a  variable  z  that  is  a  function  of  vari¬ 
ous  independent  variables 


z  =  f{xi,i:l  ^  N),  (11) 


the  uncertainty  ofz,  Az,  can  be  expressed  as  [li  1,'’)51 


/_|j;(443K^j,(490)_'| 
Vr„(4)  +  5^gg{r„(667)A 

a(/lo)  =  a„(/lo)  +  10-1  i-'e-i.aeer-o. 469/ 


(5) 


Az  = 


N 


^(Az(Ajc,) 


(12a) 


where  a^jiAo)  is  the  contribution  of  pure  water  [•]•]]. 
Then  6fop(Ao)  is  calculated  by 


^bpi^o)  —  ^{^0)  a(2o)  “  bbt^{Ao),  (6) 

To  calculate  a{A)  requires  extrapolation  of  at  Aq 
to  A.  The  relationship  between  6^^  at  two  wave¬ 
lengths  can  be  generally  expressed  as 


bbpW  =  bbp{Ao)  p(A,Aq).  (7a) 


For  easy  extrapolation  and  following  tradition 
[25,2  1],  p(A,  Aq)  in  QAA  is  modeled  as 

p(^,4)=(yy',  (7b) 


with  /;  estimated  from 


n 


-2.0 


1  -  1.2  exp 


A-rs(440)\ 

rrs(550)y 


(8) 


Here  Az(AjCj)  is  the  uncertainty  inz  that  results  from 
uncertainty  ofxj,  Axj.  This  equation  is  also  commonly 
expressed  as  [21] 


(12b) 


with  dz/dxi  a  partial  derivative  of  z  over 

In  the  QAA  process,  the  variables  rr3(A),  a(Ao), 
p{A,Aq)  (or  //),  f,  and  ^  are  involved  at  various  steps 
(see  Fig.  1).  Although  their  values  may  be  estimated 
from  the  same  rrs(A),  e.g.,  a(Ao)  or  ^/,  there  are  no  ob¬ 
servations  indicating  that  their  uncertainties  covary; 
therefore  the  uncertainties  of  these  variables  are 
considered  independent  of  each  other.  Also,  it  is 
necessary  to  keep  in  mind  that  an  lOP  product  is 
not  necessarily  affected  by  all  of  these  variables. 
For  instance,  the  value  of  p(A,Ao)  (or  //)  has  no  effect 
on  the  estimation  of6^,p(Ao),  and  the  values  of  f  and  ^ 
have  no  effect  on  the  estimation  of  a  (A).  We  recognize 
that  r^s  from  satellite  or  any  platform  is  not  error  free 
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(e.g.,  Antoine a/.  [7],Hue^a/.  [oG],  and Zibordie^ a/. 
[7,37]),  but  this  measurement-introduced  uncer¬ 
tainty  is  omitted  in  this  study  as  it  is  not  yet  clear 
about  the  pixel-specific  uncertainty  of  a  measured 
The  following  is  thus  focused  on  how  the  un¬ 
certainties  of  a(ylo),  affect  the  quality  of 

derived  lOPs. 

a.  Uncertainty  in  estimated  6^p(A):  When 
Aa(/lo)  is  known,  the  uncertainty  of  based 

on  Eqs.  (Tia)  and  (b),  is 

(13a) 

assuming  B  is  error  free.  Equation  (Iba)  indicates 
that  the  uncertainty  of  b^p{XQ)  is  proportional  to  that 
ofa(Ao). 

The  uncertainty  of  at  other  wavelengths,  how¬ 
ever,  will  be  different  because  it  depends  on  the 
accuracy  of  see  Eq.  (7a).  When  p{X,Xq)  is 

described  by  Eq.  (71)),  the  uncertainty  of  bijp{X)  then 
depends  on  the  accuracy  of  //: 


Equations  (laa),  (h’^d)),  and  (10  indicate  that,  in 
the  QAA  system,  when  both  a(Ao)  and  tj  are  error 
free  [both  Aa(Ao)  and  At]  are  zero],  the  derived  back- 
scattering  coefficient  and  absorption  at  other  wave¬ 
lengths  will  also  he  error  free.  This  suggests  that, 
in  addition  to  improving  the  qualities  of  measure¬ 
ments  and  r^s  models,  we  should  focus  on  minimizing 
Aa(Ao)  and  At]  in  future  efforts  to  improve  QAA 
performance. 

c.  Uncertainties  in  estimated  aph{X2)  and  ^^^(^2): 
Subsequently,  we  can  analytically  evaluate  the 
propagated  uncertainty  in  QAA-derived  aph{A2)  and 
^dg{^2)>  respectively.  Since  a^(A)  are  considered  con¬ 
stants,  there  are  four  properties  [see  Eq.  ( I  ())]  contri¬ 
buting  to  the  uncertainties  of  QAA-derived  ap/t(A2) 
and  adg{X2):  a(Ai),  ci{X2),  The  uncertainty 

contributions  of  f  and  ^  are  treated  independent  of 
each  other,  as  field  observations  have  not  demon¬ 
strated  that  C  covaries  with  f  in  natural  aquatic 
environments. 

The  uncertainty  contributions  of  a(>li)  and  0(^2) 
to  Aa(ig{X2)  and  Aap^{X2)  are  through  the  partial 


+  (^[5(^)  a{Ao)  -  66,„(Ao)]  In  A,,)'. 


(13b) 


In  Eq.  (]bl>),  the  first  term  on  the  right-hand  side  is 
for  At]  =  0,  while  the  second  term  is  for  Aa(Ao)  =  0- 

Note  that  Eq.  (71))  is  not  error  free  in  describing  the 
spectral  variation  of  6fcp(/l),  and  this  imperfection  is 
lumped  into  the  uncertainty  of  variable  t].  For  two 
different  wavelengths,  e.g.,  411  and  443  nm,  it  was 
found  that  /;(411,555)  is  highly  correlated  with 
/>(443,555)  [i?2  =  o.93,  NOMAD  (the  NASA  bio- 
Optical  Marine  Algorithm  Dataset  [bS]),  not  shown 
here].  It  is  thus  justified  to  use  the  combination  of 
Eq.  (71))  and  At]  for  the  estimation  of  A/;(A,Ao)  and 
then  A6^p(/l). 

b.  Uncertainty  in  estimated  a  (A):  The  uncertain¬ 
ty  of  a  (A)  resulting  from  Aa(Ao)  and  At]  is  derived  by 
combining  Eqs.  (9)  and  (13b),  omitting  uncertainties 
associated  with  quantities  A  and  B,  which  are  calcu¬ 
lated  from  measured  Trs  through  Eqs.  (1)  and  (2): 


differences  between  a(Ai)  and  a(A2)  (a(Ai)-^a(A2) 
and  ^a(A2)-a(Ai),  respectively;  see  Eq.  (10)).  There¬ 
fore,  it  is  necessary  to  know  the  algorithm-introduced 
uncertainty  of  the  partial  differences  for  analyzing 
Aarf^(A2)  and  Aap/j(A2).  Define  a  12;  as 


ai2c  =  a{X,)-Ca{X2).  (15) 

Then,  after  a  simple  mathematical  manipulation. 


^12;  —  bbp(^o) 


+  (16) 


Here  a  is 


Aa(A)  = 


A(A)B(Ao)(^ 


Aa(Ao) 


A(A)[B(Ao)a(Ao)-6fc„(Ao)](^ 


A,, 


(14) 
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which  has  no  dependence  on  both  a(/^o)  and  tj. 

Following  Eq.  ( 1 21)),  when  bfyp{^,Q)  and  tj  have  uncer¬ 
tainties,  the  uncertainty  of  ai2<;,  is 


+  \  bhp i^o) 


-CA(A2) 

Ay.)|iV 


SI}’ 


-CA(;2)(^y'ln 


SIM 


(18) 


Therefore,  with  and  for  uncertainties  of  f 
and  <J,  respectively,  Aadg{^2)  is 

(Aa,,(A2))^  = 

Qi2i;  ^ 


ciwih)  -  ci{h) 


AC 


^12C 


AC  .  (19) 


with  inverted  aph{^.2)  a^^(/^2)  ^re  much  more 

complex  when  compared  to  that  of  bfyp{X)  or  a(A), 

3.  Data  to  Quantify  Uncertainties  of  Inherent  Optical 
Properties  Derived  by  Quasi-Analytical  Algorithm 

Ideally,  error-free  field  measured  data  are  used  to  de¬ 
rive  and  quantify  uncertainties  of  algorithm-derived 
products.  When  an  error-bearing  data  set  is  used,  it 
is  difficult  to  assign  the  difference  between  algorithm 
outputs  and  measurements  to  algorithm  error  or  mea¬ 
surement  error.  Because  there  is  no  error-free  data  set 
from  field  measurements,  and  the  purpose  of  this 
study  is  to  isolate  and  evaluate  uncertainties  intro¬ 
duced  by  the  algorithm  alone,  we  simulated  a  data 
set  that  best  represented  natural  variability,  includ¬ 
ing  both  lOPs  and  r^s,  in  a  fashion  similar  to  the  Inter¬ 
national  Ocean-Colour  Coordinating  Group  (lOCCG) 
Algorithm  Working  Group  [19]. 

Generation  of  the  data  set  starts  with  the  creation 
of  spectral  lOPs,  6^  and  a,  in  particular: 

a(^)  =  +  ap/i(A)  + 

bbW  =  bf,M)  +  bf,p{X).  (23) 

Values  for  and  65^^ (>^)  are  already  known 

[.)3,39],  and  the  following  optical  models  were  used 
to  create  spectra  of  the  other  components: 


Similarly,  for  the  uncertainty  of  ap/j(A2),  we  define 
021;  as  021;  =  ^a(/l2)  -  o(Ai). 

Thus,  Aaph{^2)  Is 


(Ao,v,(a2)) 


2 

(s'-O" 
+ 


021;  +  Q!f('^l)  ~  <?Qw('^2)  .  ^ 

\ 

Q21f  +  ~  ^0-u){^2) 


,  (21) 


with 

(Aa2i;)^  =  \  ^bhpiAo) 


^(^2)(^)"-A(1i 


^bp  {-^0  ) 


^(^2) 


-A(Ai 


/ln\ '' 


In 


A,, 


I’ 


(22) 


In  Eqs.  ( i  9)  and  (li  1 ),  values  of  all  parameters  are 
given  or  derived  from  r\Q{A);  consequently  the  uncer¬ 
tainties  of  QAA-derived  cidg{^2)  ciph{^2)  t)e  es¬ 

timated  for  each  rrs(>^).  Because  api^{A2)  and  cidg{A2) 
are  next  level  products  from  r^g  inversion  [2()],  i.e., 
derivative  products  of  a  (A),  uncertainties  associated 


=  Sa(‘^40)o+^(1), 
adA^)  =  ad^(440)e-5('-44'», 

/440\ n 

fcbp(^)  =  fobp(440)(— j  .  (24) 

Here apf^i^)  is  theap/j{440)-normalized  phytoplank¬ 
ton  absorption  coefficient,  which  provides  a  spectral 
shape  forap/i(A).  Both ap/j(440),in  arange  of '^0.0056- 
0.42  m"\  and  ciL(A)  were  taken  from  the  lOCCG  da¬ 
tabase  [19].  Although  the  spectral  and  bfyp  models 
are  idealistic,  the  small  deviations  are  considered  hav¬ 
ing  negligible  effects  to  the  analyses  here. 

Because  the  estimation  of  a(Ao)  in  QAA  depends 
completely  on  the  spectral  shape  of  rrs(^)>  which  is 
primarily  determined  by  the  shape  and  magnitudes 
of  aph{A)  and  a^JA)  (they  are,  in  general,  selectively 
stronger  in  the  shorter  wavelengths)  and  the  shape  of 
bfyp{A)  (its  magnitude  is  nearly  canceled  out  in  the  Trs 
band  ratios),  we  systematically  varied  and  // 

for  each  aph{A),  but  left  65p(440)  varying  in  a  random 
way,  similar  to  the  lOCCG  data  set.  Specifically,  the 
magnitude  parameters  of  a^^,(440)  and  6^p(440)  were 
determined  as  follows: 


arf^,(440)  =piap;,(440), 

6,p(440)  =p2(ap,(440)  +a^,,(440)).  (25) 

The  value  pi,  which  describes  the  ratio  of 
a^^,(440)/ay/^(440),  was  varied  from  0.2  to  7.0  with 
a  step  of  0.2  (35  pi  values).  The  spectral  slope  of 
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adg{X),  Sy  was  varied  from  0.01  to  0.02  with  a  step  of 
0.002  (6  S  values).  The  spectral  shape  of  6^^ (A), ;/,  was 
varied  from  0  to  2.0  with  a  step  of  0.2  (11  //  values). 
The  value P2  was  used  to  generate  6^p(440)  (and  then 
bbp{X)y  when  associated  with  values  of  ;/)  and  was 
determined  by  the  following: 


P2  =  0.001  -f 


0.331a^/,(440) 
0.006  4-  a^/,("440)  ’ 


(26) 


with  31  a  random  value  between  0  and  1.  Therefore 
P2  is  in  a  range  between  0.001  and  0.3,  and  between 
0.001  and  0.15  at  the  lower  end  and  a  wider  range  for 
larger  ap;j(440)  values,  in  a  way  to  mimic  natural 
variations.  In  QAA,  because  b^pi^)  is  analytically  de¬ 
rived,  the  way  of P2  variation  does  not  matter  here,  as 
long  as  the  overall  range  is  consistent  with  natural 
environments. 

Based  on  the  above  setup,  there  are  2310  (35  x  6  x 
11)  combinations  ofaj^(/l)  and  b^pi^)  for  each  (440) 
value.  Since  there  are  20  different  ap/j(440)  values  in 
the  lOCCG  database,  we  created  a  data  set  with 
46,200  combinations  of  cLdg{^)y  and  bbp{^)y 

and  then  a{^)  and  6t(/l).  Note  that  is  not  the 

same  for  the  same  ap;j(440)  value  in  the  lOCCG  data¬ 
base;  thus  the  spectral  variability  of  aph{A)  within 
each  group  is  maintained  here. 

After  a  (/I)  and  are  known,  to  be  more  efficient, 

rrs(A)  is  modeled  using  Eq.  (2)  instead  of  using  Hydro¬ 
light  [10].  This  is  also  because  we  are  assuming  a  per¬ 
fect  analytical  model  between  the  apparent  optical 
properties  (AGP:  rrs(yl))  and  the  lOPs  (6^  and  a)  and 
focusing  on  uncertainties  of  a(>^o)  ;/  (and  f  and 

(J)  and  their  propagations  in  the  QAA  system.  Indeed, 
there  is  a  <10%  error  associated  with  Eq.  (2)  [29,1 1], 
which  has  a  bigger  effect  on  the  retrieval  of  b^p  when 
bfyp  is  extremely  small,  e.g.,  for  waters  in  the  oceanic 
gyres  [12].  However,  errors  associated  with  the 
model  can  be  corrected  later  when  a  more  accurate 
model  is  adopted  (Lee  et  a/,  in  preparation),  and  such 
a  correction  will  not  affect  the  general  conclusions 


about  the  uncertainties  of  a(2o),  f]y  f  >  and  ^  and  their 
propagations  to  other  lOPs. 

The  wavelengths  used  in  this  study  (410,  440,  490, 
550,  and  670  nm)  do  not  exactly  match  those  from  the 
current  operational  satellite  sensors  such  as  Sea- 
WiFS/MODIS/MERIS.  The  spectral  differences  of  the 
corresponding  bands,  however,  are  so  small  that  re¬ 
sults  and  conclusions  achieved  here  are  applicable  to 
those  ocean  color  systems.  Figure  2(a  >  presents  values 
and  ranges  of  the  simulated  a(440)  and  a(550),  and 
Fig.  2(bi  illustrates  values  and  ranges  of  rrs(440) 
and  ri.s(550).  The  figures  also  show  equivalent  data 
from  field  measurements  (NOMAD)  [38].  Both  charts 
demonstrate  a  wide  range  of  values  at  440  nm  for  a 
given  value  at  550  nm,  as  expected  in  natural  aquatic 
environments  (simulated  absorption  actually  has  a 
wider  range  than  current  NOMAD  values).  These 
comparisons  indicate  that  the  simulated  data  set  is 
consistent,  in  both  magnitude  and  variability,  with 
data  from  field  measurements.  A  unique  feature  of 
the  simulated  data  set,  however,  is  that  it  is  free  of 
measurement  errors,  and  thus  well  suited  to  identify 
and  analyze  algorithm  introduced  uncertainties. 

4.  Results  and  Discussion 

A.  Uncertainty  in  the  Estimated  a(4) 

As  shown  by  Eq.  (11),  one  key  component  in  the  QAA 
system  is  the  estimation  ofa(/lo),  because  it  is  one  of 
the  properties  that  initiates  the  entire  sequential  pro¬ 
cess.  Figure  3(a)  illustrates  derived  a(550)  [from 
Eq.  (5),  550  nm  is  considered  Aq  here;  also  note  that 
sensor  specific  coefficients  are  required  to  obtain  con¬ 
sistent  results  across  different  sensors]  versus  known 
a(550).  Because  Eq.  (5)  is  not  perfect,  there  is  an  error 
(big  or  small)  associated  with  each  inverted  a (550).  To 
highlight  the  range  of  errors  between  known  and 
inverted  a(550).  Fig.  3(1))  shows  the  distribution  of 
the  absolute  percentage  error  (ape):  |a(550)i^o^j^- 
a(550)Qy^|/a(550)i^Q^j^.  About  57%  of  the  data  have 
an  ape  within  '-^13%,  and  about  70%  of  the  data  have 
an  ape  within  20%,  The  overall  average  of  ape  is 
15.6%. 


(a) 

LvLj'.IL'a  ' 

I 

Siniula<K>n  (462O0) 

•  NOMAD  (91 S) 

00  02  0.4  06  0.8  1.0 

^/(550)  [m  '] 


/- ,(550)  [sr  ^] 


Fig.  2.  Examples  of  simulated  optical  properties  at  440  and  555  nm  (a)  for  absorption  coefficients  and  (b)  for  remote-sensing  reflectance 
(subsurface).  For  comparison,  corresponding  values  from  the  NOMAD  data  set  [.IS]  are  also  plotted.  For  NOMAD,  there  are  915  pairs  of 
absorption  coefficients,  while  there  are  984  pairs  of  reflectance  coefficients,  but  there  are  46,200  pairs  of  values  for  the  simulated  data  sets. 
The  NOMAD  values  were  converted  from  above-surface  remote-sensing  reflectance  with  =  Rf^/{0.52  +  l.7R^). 
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Fig.  3.  (a)  QAA-dorived  a(550)  compared  with  known  a (550)  of  the  simulated  data  set.  Statistics  of  data  in  log  scale  are  shown  in  the 

figure,  (b).  Histogram  of  absolute  percentage  error  in  linear  scale,  |a(550)j^,,o,^n  -  a(550)Qy^  | /a  (550)  "The  average  of  this  error  is  15.6%. 


To  evaluate  the  error  for  different  a  (550)  values,  we 
divided  the  data  set  into  30  groups  that  correspond  to 
a  narrow  range  ofa(550)  (within  ±2%  of  a  designated 
a(550)).  The  designated  a(550)  values  are  0.058, 
0.063, 0.069,  ...,0.374, 0.389,  to  0.405  with  an  in¬ 
crease  rate  of  ^8%  at  the  lower  end  and  ^4%  at  the 
higher  end.  As  remotely  the  “true”  a (550)  is  unknown, 
the  designation  ofa(550)  is  based  on  QAA-derived  va¬ 
lues.  Note  that  the  a (440)  is  in  a  range  of  ^0.013- 
0.085  for  the  group  with  a(550)  =  0.058 
(slightly  higher  than  the  absorption  coefficient  of  pure 
water  at  this  wavelength  [o3l,  0.0565  m"^),  and  is  in  a 
range  of  ^0. 9-2.8 m"^  for  the  group  with  a(550)  = 
0.405  m”h  Such  ranges  indicate  that  the  a(550) 
range,  ^0.058-0.4  covers  all  oceanic  and  most 
coastal  waters.  We  limited  the  analyses  to  a(550)QAA 
<  ^0.4  m“^  also  because  that  QAA-derived  a(550)  gra¬ 
dually  plateaus  at  ^^0.5  m”\  which  suggests  that  the 
current  version  of  QAA  contains  higher  uncertainties 
when  QAA-derived  a(550)  >  0.5  m”\  i.e.,  for  extre¬ 
mely  turbid  waters.  For  such  cases,  a  reference  wave¬ 
length  longer  than  550  nm  is  needed  (e.g.,  Lee  et  ai 
123]  and  Doron  et  aL  [11]). 

For  each  point,  the  absolute  difference  (f),  or 
error,  between  QAA-derived  and  known  a (550)  is 
calculated: 


f  ^  Ia(550)knowi,  -  a(550)QAAl,  (27) 


along  with  its  relative  difference,  i.e.,  ^*/a(550)Q^. 
As  examples.  Fig.  1  shows  the  distributions  of  e/a 
(550)q^  for  six  groups  with  different  a(550)Q^  val¬ 
ues.  For  this  simulated  data  set  with  the  current  ver¬ 
sion  of  QAA,  the  relative  differences,  f/a(550)Q^, 
are  centered  at  the  lower  end  of  the  values.  In  parti¬ 
cular,  for  waters  with  low  a(550)Q^(<  0.076 m"^), 
70%  or  more  of  the  true  a (550)  are  within  ±10%  of 
a(550)Q^;  when  a(550)Q^  becomes  larger,  e.g., 
0.2  m"^  or  more,  some  of  the  true  a(550)  could  be 
off  by  ±60%  of  a(550)Q^.  These  results  demonstrate 
that  the  errors  or  uncertainties  in  ocean  color  remote 
sensing  are  not  uniform,  as  recently  found  by  Moore 


et  aL  [9]  when  empirically  retrieved  chlorophyll-a 
concentrations  were  evaluated. 

The  uncertainty  of  QAA-derived  a(550),  Aa(550), 
is  defined  as  the  average  of  e  for  each  group  because 
it  approximates  the  65th  percentile  of  the  r  distribu¬ 
tion;  see  Fig.  By  this  definition,  there  is  a  65% 
likelihood  that  the  “true”  a (550)  will  be  within 
a(550)Q^  ±  Aa(550).  Figure  5  shows  the  relative 
uncertainties  (Aa(550))/a(550)Q^,  with  u(550)q^ 
the  average  a(550)Q^^  of  each  group)  for  the  various 
a (550)  groups.  Fora(550)Q^  <  0.1  m“^  which  covers 
^^95%  of  global  waters  (Bryan  Franz,  personal  com¬ 
munication),  Aa(550)/a(550)Q^  is  in  general  <16%. 
Since  a  (550)  is  dominated  by  tne  contribution  of  pure 
water,  the  uncertainty  associated  with  QAA  derived 
a(550)  is  limited  for  most  oceanic  waters.  On  the 
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Fig.  4.  Histograms  of  |a(550)^o^,,  -a(550)QAAl/a(550)t^AA  ^ 

few  designated  a(550)Qy^  values.  The  first  value  in  a  figure  repre¬ 
sents  a(550)Q^,  and  the  second  value  represents  the  number  of 
points  for  that  a(550)Q^  group  (see  text  for  details). 
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group.  The  solid  dots  represent  an  empirical  fit  [Eq.  (2.S)  divided 
by  a(550)Q^l  for  easy  evaluation  of  Aa(550). 
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other  hand,  uncertainties  are  higher  [larger  Aa 
(550)]  for  higher  a  (550)  values,  where  marine  optical 
and  biogeochemical  properties  are  more  complex. 

To  derive  Aa(550)  for  any  given  a(550)Qy^,  a  best- 
fit  empirical  relationship  between  Aa(550)  and 
a(550)Q^  was  derived  (dotted  line  in  Fig.  5); 

Aa(550)  =  0.35(1  -  2.4exp(-16.0a(550)Q^)) 

xa(550)QAA-  (28) 

Because  a (550)  in  our  simulated  data  set  has  a 
wider  range  and  higher  variability  than  in  the  NO¬ 
MAD  data  set  (see  Fig.  2),  Aa(550)  of  real  environ¬ 
ments  could  be  smaller  than  that  represented  by 
Eq.  (2S),  assuming  high-quality  r^g. 

In  the  above  analysis  the  number  of  points  of  each 
group  spans  from  2812  fa(550)QAA  =  0.058  m"^]  to  401 
[a(550)Q^  =  0.405  m"^].  To  test  the  effect  of  the  num¬ 
ber  of  points  on  the  evaluation  of  Aa  (550),  we  changed 
the  lower-upper  bounds  of  each  group  from  within 


±2%  to  within  ±5%  of  the  designated  a  (550),  result¬ 
ing  in  nearly  double  the  number  of  points  for  each 
group  [4873  for  a(550)Q^  =  0.058 m“\  and  962  for 
a(550)Qy^  =  0.405  m"^].  We  found  negligible  effects 
to  the  statistical  results  shown  in  Fig.  5,  which  pro¬ 
vides  confidence  that  Aa(550)  values  shown  in  Fig.  5 
are  not  sensitive  to  a  (550)  grouping  so  long  as  narrow 
bounds  are  maintained. 

B.  Uncertainty  in  the  Estimated  /^^^(Ao) 

With  Aa(550)  known,  the  uncertainty  of  5^p(550)  can 
be  easily  calculated  [see  Eq.  (13a)].  Figure  G{a)  pre¬ 
sents  the  relative  uncertainties  of6ip(550)  [A5^p(550) 
/6^p(550)qy5^]  of  the  simulated  data  set,  which  are 
generally  in  a  range  of  ^-^2-40%,  but  can  approach  in¬ 
finity  when  derived  6ip(550)  nears  zero.  Although 
Aa(550)  is  fixed  for  a  given  a (550 )Qy5^,  the  relative  un¬ 
certainty  of  b^p  (550)  differs  for  the  same  QAA-derived 
6^p(550)  values.  This  can  be  explained,  combining 
Eqs.  (())  and  (13a),  by 


A6^p(550)  B(550)Aa(550) 

^>6p(550)qaa  “  B(550)a(550)QAA  "  fc6».(550) ' 


Different  combinations  ofB(550)  anda(550)Q^  can 
result  in  the  same  6i^(550).  Different  a(550)Q^  val¬ 
ues,  however,  come  with  different  Aa  (550)  values,  and 
B(550)  does  not  cancel  out  in  Eq.  (29)  until  6^^, 
(550)  is  negligible,  which  therefore  results  in  different 
A6i^(550). 

To  illustrate  how  QAA  performs  for  oceanic  waters. 
Figure  Gib)  shows  the  relative  uncertainties  of 
6^p(550)  for  data  with  a(550)  ^  0.065  m"^  and  indi¬ 
cates  that  generally  the  relative  uncertainty  is  less 
than  or  around  ±10%  for  5^p(550)Qy5^  in  the  range 
of  0.0003-0.003  m"^.  Such  a  result  suggests  that  the¬ 
oretically  highly  reliable  6^p(550)  can  be  retrieved 
from  spectral  remote  sensing  reflectance  for  such 
waters  by  QAA.  For  realistic  systems,  however,  the 
accuracy  of  6^p(550)  also  depends  on  the  accuracy 
ofB(550)  and  the  absorption  [33,43]  and  backscatter- 
ing  coefficients  [1  1,45]  of  pure  water. 
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C.  Uncertainty  in  the  Estimated  «(/l) 

The  calculation  of  Aa(yl)  requires  knowledge  of  the 
uncertainty  of  which  is  represented  by 

A//  [Eq.  (M)].  Unfortunately,  the  spatial  variability 
of  A;;  is  unknown,  although  global  maps  of  //  have 
been  produced  [2;^  !^].  Also,  there  is  not  enough 
high-quality  measured  data  to  evaluate  a  confidence 
range  of  //  for  various  regions,  although  the  expres¬ 
sion  [Eq.  ((S)]  generally  captures  the  trend  of  higher 
values  in  oceanic  waters  and  lower  values  in  coastal 
waters  [47].  Unlike  the  total  absorption  coefficients 
that  have  a  general  spectral  dependence,  governed 
by  the  spectral  shapes  of  the  multiple  components, 
}]  is  a  parameter  used  to  determine  the  spectral  de¬ 
pendence  of  b^p.  Therefore,  unlike  Aa(4o),  A//  could 
not  be  evaluated  with  a  simulated  data  set.  Instead, 
based  on  a  comparison  (not  presented  here)  between 
Trs  estimated  tj  [Eq.  (S)]  and  calculated  //  from  mea¬ 
sured  biyp{A)  of  the  NOMAD  data  set,  a  A^  =  0.5, 
which  represents  the  68th  percentile  of  the  absolute 
tj  difference,  was  assigned  universally.  This  assigned 
value  represents  a  high  variability;  i.e.,  if  the  true  // 
value  is  1.0,  A;/  =  0.5  suggests  an  //  range  of  0.5-1. 5 
[/)(440. 550),  range  of  '^1.12-1.40]  for  68%  of  the  sam¬ 
ples,  which  cover  nearly  all  waters. 

To  show  the  different  contributions  of  Aa(4o)  and 
Atj  to  Aa(4),  Eq.  (1  I)  is  applied  to  the  simulated  data 
set.  As  an  example,  Aa(440)  is  calculated  with 
Aa(550)  and  A//  separately,  and  then  calculated  with 
both.  Figure  7  presents  their  relative  contributions 
and  the  combined  effects,  with  the  Y  axis  for  Aa (440)/ 
a(440)Q^,  the  relative  uncertainty.  Clearly,  1)  there 
are  different  Aa(440)  for  the  same  a(440)Q^  (see 
Fig.  7),  as  its  inversion  is  also  associated  with  particle 
scattering  properties,  and  2)  Aa(550)  and  A;/  con¬ 
tribute  differently  to  Aa(440).  For  relatively  clear 
waters  [a(550)  <  0.08m”^],  because  a(550)Q^  is 
quite  accurate,  the  contribution  of  Aa(550)  to 
Aa(440)  is  limited  (the  blue  dots  in  Fig.  7).  The  con¬ 
tribution  of  A;/  to  Aa(440),  on  the  other  hand,  could  be 


Fig.  7.  (Color  online)  Relationship  between  relative  uncertainty 
of  a(440)[Aa(440)/a(440)Qy^,  red]  and  QAA-derived  a(440).  Also 
shown  are  contributions  of  Aa(550)  (blue)  and  At]  (green), 
respectively,  to  the  relative  uncertainty. 


higher  than  that  from  Aa(550),  as  its  effect  also  de¬ 
pends  on  the  value  of  B(550),  which  is  directly  related 
to  the  value  of  ri.s(550)  and  the  model  used.  For 
coastal  turbid  waters,  however,  the  maximum  contri¬ 
bution  from  A//  is  ^11%,  which  approximates  the 
square  root  of 550/440.  But  Aa(550)  in  general  contri¬ 
butes  more  to  Aa(440)  for  complex  waters,  simply  be¬ 
cause  Aa(550)  itself  gets  larger.  Note  that  for  a 
specific  case,  because  the  overestimation  or  underes¬ 
timation  ofa(4o)  or  //  are  not  necessarily  in  phase  for  a 
measured  rrs(4),  the  effects  of  errors  in  a(4o)  and  //  on 
^bp  ('^)  a(4)  are  not  necessarily  compounded.  If  one 

is  overestimated  while  the  other  is  underestimated, 
the  actual  error  in  the  derived  b]jp(A)  and  a(4)  could 
be  minimal. 

The  above  results  suggest  that  for  waters  with 
a(440)Q^  <  0.05  m"^  (most  oceanic  waters,  Bryan 
Franz,  personal  communication),  the  “true”  a(440) 
value  is  most  likely  within  d=9%  of  a(440)Q^,  and 
could  be  within  ±2%  of  ct(440)QAA»  ^'rs  measure¬ 
ments  and  models  are  highly  reliable.  This  result 
and  the  result  about  6^p(550)  indicate  that  since  r^s 
is  a  cumulative  measure  of  the  upper  tens  of  meters 
of  such  waters  [  IH,  inversion  is  an  efficient  and  re¬ 
liable  sampling  method  for  measurement  of  optical 
properties  of  many  oceanic  waters.  Note  that  tradi¬ 
tional  water  sampling  approaches  require  either  long 
path-length  tubes  (e.g.,  ac-9,  Wetlabs,  Inc.)  or  a  large 
volume  of  water  samples  [19,50]  to  obtain  reliable  in 
situ  measurements  for  such  waters. 

Separately,  because  the  diffuse  attenuation  coeffi¬ 
cient  of  downwelling  irradiance  (K^)  is  a  simple  func¬ 
tion  of  a  and  6^  [5 1 ,52],  the  above  analysis  can  easily  be 
extended  to  evaluate  the  uncertainties  associated 
with  the  semianalytically  derived  Kfi  products  [5 1 ,52] . 

D.  Uncertainty  in  the  Estimated  ^^^.(440)  and  a^;,(440) 

The  uncertainties  of  a^(440)  and  ap/j(440)  are  asso¬ 
ciated  with  both  the  uncertainties  of  the  QAA-derived 
a  (4)  and  the  uncertainties  of  parameters  ^  and  <J.  The 
parameter  (  varies  in  a  relatively  small  range  (0.7- 
1.0)  based  on  field  measurements  [52],  but  ^  could 
have  much  wider  variations.  When  fitting  measured 
a^,(4)  with  an  exponential  function,  the  spectral  slope 
(S)  is  generally  in  a  range  of  0.01-0.02  nm”^  [49,54- 
5()],  which  suggests  a  range  of  ^1.35-1.82  for 
Because  it  is  not  clear  yet  how  to  best  assign  values 
for  both  (  and  ^  for  an  individual  pixel,  we  chose  to 
assign  a  constant  value  of  0.85  for  ( (the  mode  of  NO¬ 
MAD,  not  shown  here)  and  1.568  for  which  is  equiva¬ 
lent  to  S  =  0.015  nm"^,  or  a  power  value  of  6.8  when 
adg{A)  is  fitted  with  a  power-law  function  of  wave¬ 
length  [5(i].  Further,  we  set  a  universal  uncertainty 
of  A(  as  0.1  and  A^  as  0.14,  which  indicates  a  (  range 
of  0.75-0.95  and  a  ^  range  of  1.43-1.71  (equivalent  S 
range  of  0.012-0.018  nm"U;  both  cover  a  majority  of 
ocean  waters.  Note  that  these  values  and  ranges  could 
be  refined  later  after  we  acquire  additional  knowledge 
about  their  ways  of  variations. 

With  the  above  values  for  C,  <?,  Af,  and  Aa(550) 
and  At]  for  each  rj^{A),  values  of  Aa^(440)  and  Auph 
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Fig.  8.  (Color  online)  Uncertainty  of  QAA-derived  aj^(440).  Also 
shown  are  the  contributions  of  Aa,  A^,  and  A^,  respectively. 

(440)  were  calculated  and  are  presented  in  Figs,  s 
and  a The  separate  contributions  of  Af,  Af , 
and  Aa  are  also  presented  in  order  to  demonstrate 
their  relative  importance  in  Aa^^(440)  and 
Aap/j(440),  respectively.  We  do  not  present  percen- 
tage-like  error  for  arf^(440),  or  ap;j(440),  as  shown 
in  Figs.  T)-?,  because  when  a (410)  and  a(440)  en¬ 
counter  larger  errors,  it  can  numerically  result  in  ne¬ 
gative  or  close  to  zero  values  (^8%  of  this  data  set),  or 
very  large  values,  for  a^^(440)  [or  (440)1  (com¬ 
monly  referred  to  as  nonvalid  retrieval  in  semiana- 
lytical  algorithms  [‘J(),‘22,271).  Dividing  Aa^^(440) 
(or  Aap/j (440)1  to  such  QAA-derived  a^^(440)  [or 
Up/j  (440)1  values  could  result  in  a  misleading  per¬ 
spective  of  relative  uncertainty.  The  Aa^^(440)  [or 
Aup/j  (440)1  values,  along  with  their  corresponding 
a^^,(440)  [or  Up/j  (440)1  values,  provide  adequate  infor¬ 
mation  about  the  reliability  of  a^^(440)  [or  ap/j(440)] 
derived  from  a  spectral 

A  few  important  findings  can  be  drawn  from  the 
results  shown  Figs.  <S  and  0(a): 

a.  Not  surprisingly,  Aa^^(440)  and  Aap;j(440) 
vary  a  lot.  Their  values  and  distributions  are  differ¬ 
ent  from  each  other,  and  different  for  the  same  QAA- 
derived  af/^(440)  and  ap;j(440)  values,  respectively. 


i.e.,  nonuniform  uncertainties.  And,  for  the  same 
QAA-derived  a^^,(440)  and  ap;i(440)  values,  Aa^^ 
(440)  is  apparently  smaller  than  Aap;j(440).  For  in¬ 
stance,  when  QAA-a^^(440)  is  0.1  m"^  ^a^^{440) 
could  be  as  low  as  0.02  m"\  or  as  high  as  0.09  m“^ 
When  QAA-ap;j(440)  is  0.1  m“\  however,  Aap;i(440) 
could  be  in  a  range  of  0.02-0.7  m“^  Such  a  result  sug¬ 
gests  that  fundamentally  we  could  do  better  in  anal¬ 
ytically  retrieving  a^g  than  retrieving  apf^.  The 
results  also  highlight  the  importance  and  necessity 
of  obtaining  quality  measurements  for  each  inverted 
a^g  and  aph  from 

D.  Between  Aa^^(440)  and  Aap;j(440),  it  is  inter¬ 
esting  that  Aa^«(440)  generally  increases  with 
QAA-derived  arf^(440),  but  not  QAA-derived  Aup/^ 
(440).  This  is  primarily  because  the  inversion  of 
a^g{440)  is  proportional  to  a(Ai)-Ca(A2),  while  inver¬ 
sion  of  ap;j(440)  is  proportional  to  [see 

Eq.  (10)1.  When  absorption  in  the  blue  domain  is  ac¬ 
tually  composed  with  high  a^ig  that  is  also  associated 
with  a  large  ^  value,  say  ^l.§  (equivalent  S  value  of 
0.019  nm”^)  and  thena(Ai)  ^  1. 8a (^2),  it  may  result  in 
very  small  or  negative  ap/j(440)  when  a  value  of 
1.568  (equivalent  S  of  0.015  nm"^)  is  used  in  QAA 
(or  other  model-based  analytical  inversion  methods). 
For  such  a  case,  although  the  QAA-derived  ap;i(440) 
could  be  near  zero  or  negative,  Aap;i(440)  could  be 
very  large,  from  a  large  Aa. 

c.  Obviously  Af,  A^,  and  Aa  all  contribute  to 
Aa^^(440)  and  Aap/j(440).  It  is  worth  noting,  how¬ 
ever,  that  A<J  and  Aa  generally  contribute  more  than 
A(^.  In  most  cases,  Aa  values  contribute  more  than 
A^.  This  is  because,  since  ^  is  much  larger  than 

a  small  change  in  a  (A)  could  more  strongly  affect 
ap;j(440).  Sometimes  A^  could  contribute  more  than 
Aa  to  Aap;j(440)  ,  again  when  a^^  dominates.  These 
results  suggest  that,  as  a  first  order  estimation,  it  is 
justified  to  use  universal  default  f  and  ^  values  for 
operational  processing  ocean  color  images  of  global 
oceans.  For  further  improvement  in  analytically  de¬ 
riving  Up/j  from  Trs,  it  is  more  important  to  get  a  better 
estimation  of  ^  (or  S)  rather  than  f. 

d.  Figure  9{  a )  also  highlights  the  potentially  high 
uncertainties  associated  with  ap/j(440)  when  it  is 
algebraically  derived  from  the  olue  portion  of  an 


Fig.  9.  (Color  online)  Uncertainty  of  QAA-derived  ap/,(440)  (a)  for  the  entire  simulated  data  set,  with  contributions  of  Aa,  AC,  and  A^, 
respectively;  (b)  for  data  with  0.06  5  a^^(440)/ap/,(440)  ^  2  and  1.43  ^  (4 10) (440)  <  1.72. 
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Trs  spectrum,  especially  when  water's  absorption 
coefficient  is  dominated  by  the  contribution  from 
gelbstoff  and/or  detritus.  In  order  to  mitigate  the  un¬ 
certainties  associated  with  analytical  ap^  inversion 
for  such  environments,  other  approaches  may  be 
necessary.  These  include  utilization  of  fluorescence 
line  height  [12,57],  or  using  information  at  other 
bands,  by  incorporating  spectral  models 

[22,30,5:)]  to  add  spectral  constraints.  Each  spectral 
model  will  introduce  its  own  uncertainties  because  it 
is  still  not  certain  about  which  spectral  model  for 
should  be  used  for  a  specific  water  body 
e.  On  the  other  hand,  if  we  restrict  the  analysis  to 
data  with  0.6  <  a(f^,(440)/ap/j(440)  <  2.0  and  1.43  < 
^<1.72  (equivalent  to  0.012  <  S  <  0.018  nm"^)  i.e., 
waters  usually  encountered  in  the  open  ocean  [ocS], 
the  uncertainties  of  ap/j(440)  are  much  narrower 
[see  Fig.  9<l)i].  Such  results  support  the  strategy  of 
using  rrs('l)  in  the  blue  domain  to  analytically  derive 
ap/j(440)  for  most  oceanic  waters. 

5.  Summary  and  Conclusions 

Based  on  the  algebraic  and  stepwise  nature  of  QAA 
and  theory  of  error  propagation,  we  developed  anal¬ 
ytical  expressions  [Eqs.  (13a),  (l  ib),  (11),  (19),  (21), 
and  (2S)]  for  uncertainties  of  QAA-derived  lOPs. 
These  expressions  provide  the  basis  to  evaluate 
and  understand  the  uncertainties  associated  with 
analytically  derived  lOPs  from  remote  sensing  re¬ 
flectance.  Applying  the  expressions  to  a  simulated 
wide-range  data  set,  we  evaluated  and  illustrated 
the  uncertainties  and  their  propagations  associated 
with  the  QAA  procedure  and  products.  Specifically, 
we  evaluated  the  uncertainties  involved  with  the 
QAA-derived  a(/lo),  ai'^d  further  the  propagation  of 
Aa(Ao)  and  A/;  (when  is  expressed  as  a  power- 

law  function)  to  other  analytically  derived  lOPs.  We 
found  that  the  effects  of  these  two  on  a  (A)  are  differ¬ 
ent  with  A//,  potentially  having  a  larger  effect  than 
Aa(Ao)  for  waters  with  a  low  absorption  coefficient 
but  high  scattering  coefficients,  while  Aa(Ao)  has  a 
larger  effect  for  coastal  waters.  In  general,  with 
the  current  version  of  QAA  and  assuming  that 
measurements  and  models  are  error  free,  the  most 
likely  uncertainty  of  QAA-derived  a (440)  is  under 
±13%  for  a(440)Q^  <  0.1  m"^  but  can  be  as  much  as 
±37%  for  a(440)Q^  approaching 0.5  m"^.  However, 
such  uncertainties  are  not  necessarily  the  same  for 
waters  with  the  same  a(A)Q^  values,  since  the  uncer¬ 
tainty  of  QAA-derived  a  (A)  also  depends  upon  the 
value  of  and  higher  could  cause  larger  uncer¬ 
tainty  (resulting  from  error  in  extrapolation).  For 
oceanic  waters  where  particle  backscattering  is  low, 
such  as  oceanic  gyres  [  12],  an  uncertainty  of  ±5-10% 
in  QAA-derived  a (440)  is  expected. 

Uncertainties  of  and  ^^e  much  more  com¬ 
plex  from  ocean  color  inversion,  as  their  accuracies 
depend  not  only  on  the  accuracy  of  a  (A)  and  their  re¬ 
lative  contributions  to  a  (A),  but  also  on  their  spectral 
dependencies.  There  could  be  significantly  high 
uncertainties  associated  with  algebraically  derived 


a^g  or  Up^y  if  the  total  absorption  has  a  high  uncer¬ 
tainty  and/or  one  of  the  two  components  dominates 
the  other.  However,  if  a(f^(440)  and  ap;j(440)  are  com¬ 
parable  to  each  other,  e.g.,  the  ratio  of  one  to  the 
other  is  between  0.5  and  2,  it  is  well  justified  that 
both  properties  could  be  well  retrieved,  anal5d:ically, 
from  the  blue  portion  of  rrs(A).  But,  generally,  for  the 
same  a^g{AAQi)  and  aph{AA{))  values,  it  is  found  that 
uncertainties  associated  with  ac^(440)  are  smaller 
than  that  of  aph{AAQ). 

For  accurate  derivation  o^a^g  and  the  key  com¬ 
ponents  are  the  accurate  total  absorption  coefficients 
in  the  blue/green  domain  and  the  spectral  ratio  ofa^, 
of  the  relevant  wavelengths.  For  the  derivation  of  Upfi 
in  complex  waters  where  is  not  the  primary  con¬ 
tributor  to  the  total  absorption,  methods/approaches 
other  than  simple  algebraic  inversions  are  required 
to  improve  the  quality  of  ap^  retrieval. 

These  quality  measures  can  be  extended  to  biogeo¬ 
chemical  products,  such  as  chlorophyll  concentra¬ 
tion),  when  the  uncertainties  between  the  optical 
properties  and  biogeochemical  products  are  known 
per  pixel.  Although  the  evaluation  results,  especially 
for  Aa(550),  might  be  dependent  on  the  spectral  mod¬ 
els  as  well  as  the  ranges  and  distributions  of  data 
used,  the  ranges  of  the  parameters,  S  and  //  in  parti¬ 
cular,  of  the  simulated  data  set  are  quite  inclusive 
based  on  our  current  understanding,  and  therefore 
the  results  and  conclusions  are  general  and  applic¬ 
able  to  global  oceans.  Modifications  and  refinements 
are  possible  and  expected  when  better  parameteriza- 
tions  regarding  a(Ao),  /^(A,  Aq),  and  and  for  various 
regions  are  developed. 

This  study  shows  that  we  can  not  only  derive  lOPs 
from  the  various  ocean  color  missions,  but  also  have 
necessary  components  for  quantitatively  measuring 
the  quality  of  such  lOP  products  for  each  pixel.  The 
eventual  generation  of  lOP  quality  maps,  however, 
will  depend  on  our  understanding  and  quantification 
of  the  uncertainties  of  r^s  derived  from  various  plat¬ 
forms  over  various  waters  [7,36,3)7],  as  well  as  uncer¬ 
tainties  associated  with  models  [11,59,60].  When 
the  uncertainties  of  these  variables  are  available  pix¬ 
el  wise,  they  can  be  incorporated  following  the  theory 
of  error  propagation  for  the  generation  of  lOP  quality 
maps  from  satellite  ocean  color  remote  sensing. 
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